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PREFACE 


The primary objective of this study was to demonstrate the image quality 
achievable with current optical multiplexing techniques. If acceptable, these 
techniques will significantly reduce film consumption in the distribution of dis- 
tortion-free ERTS imagery. The study was conducted over a nine-month period 
and was restricted to 35 millimeter image formats and up to twelve channel 
multiplex recordings. The results indicate that seven carrier-modulated 
images can be overlayed on a single film frame and retrieved substantially free 
from crosstalk and dimensional distortion. Resolution was maintained at 
10 line-pairs per millimeter, and the gray scale was reproduced essentially 
linear over a 100 to 1 brightness dynamic range. Based on these results, we 
recommend that the techniques be demonstrated with seven spectral zones of 
ERTS imagery in spatial registration at full size format. 
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SECTION 1 


INTRODUCTION 


SCOPE 

The resiilts of a nine-month study to evaluate multiple image overlay methods as 
applied to ERTS (Earth Resources Technology Satellite) imagery are presented in this 
technical report. Multiple imagery is a photo-optical process for storing and retriev- 
ii^ many overlapped exposures within a single film frame . Experimental camera and 
projector setups, photographic processing techniques, retrieved image results and 
system performance characteristics, as well as the theoretical models utilized for 
this study are reported. Preliminary design requirements for a full-format multiple 
image system are also presented. 

PURPOSE 

The purpose of this study was to determine whether multiple image methods that 
1 2 

had been published and patented by the contractor offered a practical, acceptable, 
and economical solution to the problem of ERTS imagery storage and distribution. 

To be practical, the methods must lend themselves to implementation within state- 
of-the-art copying techniques and materials. Laboratory experiments that can be 
conducted only by highly trained technicians under elaborately controlled conditions 
would not be considered practical. 

To be acceptable, the multiplexing methods must 3 deld an increase in storage den- 
sity without degrading imagery. The outstanding degradations to be avoided are loss 
in resolution, dimensional distortion, brightness distortion, and crosstalk. 

To be economical, the methods must exhibit a potential for cost savings as com- 
pared to alternate solutions. Reduction in film consumption must not be outweighed 
by additional equipment requirements. 

OBJECTIVES 

The primary objective of the program was to demonstrate, by means of experimen- 
tal results, the image quality that can be obtained with current multiplexing methods. 
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Using black-and-white 35 mm transparencies as inputs, we were to investigate a 
minimum of six and maximum of twelve image channels. System resolution had to be 
maintained at 10 line-pairs/mm or greater. Since the input targets for this primary 
objective were various scenes as opposed to spectral zones of the same scene, regis- 
tration of one channel to another did not apply. 

However, a secondary objective was to demonstrate how three different spectral 
zones of a single scene could be printed out in precise registration onto black-and- 
white film and subsequently retrieved as a tricolor image with no loss in registration. 
Here, the assiunption was that the spectral zones were stored with temporal regis- 
tration in an electronic buffer as opposed to being stored as spatially separated trans- 
parencies. To achieve this objective in its broadest sense would have required not 
only spectral imagery signals, but also a scanner-modulator system, neither of which 
was readily available within the scope of this study. The demonstration was therefore 
to be conducted by synthetically generating the three zone signals by means of a com- 
puter technique and printing these signals out with a step-and-repeat photofacsimile 
device that simulated a scanner-modulator system. 

DISCUSSION OF THE PROBLEM 

The ERTS program will ultimately require extremely large amounts of photographic 
materials to disseminate imagery to various users. For a substantial number of users, 
the mensurational accuracy of the imagery is of vital concern. Thus, once a master 
print has been produced by ground processing equipment, no distortions should be 
introduced by subsequent duplication procedures. 

Ideally, if all duplication were performed by means of contact printing on stable 
photographic glass plates, a minimum of distortion would be introduced. The eco- 
nomics of that approach are clearly prohibitive for mass distribution and probably 

i 

only feasible for a small number of prints to a limited number of users. Ih fact, even 
if one accepted the more practical alternative of a thick Estar base duplicating mate- 
rial, the film budget would still be quite substantial. The question is, Can large 

amounts of distortionless imagery be disseminated in a more economical way? 

% 

These inputs were used because the contractor already had encoding cameras and 
decoding projector hardware of that format. ' Additional cost would have been incurred 
to increase the image size. 
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One way to use film more economically is to utilize its full storage capacity. 
Microfilming is a prime example of efficient and economical data storage. Since we 
assume that the ERTS master prints will have an upper resolution limit of about 
10 line-pairs/mm, and films with resolving power far in excess of 500 line-pairs/mm 
are common, film consumption could be drastically reduced by duplicating at a 50:1 
reduction ratio. Lenses are readily available with the required resolution and format 
capacity. Unfortunately, only a symmetric lens used to image 1:1 at fixed conjugates 
can be made to operate substantially free of distortion. 

Thus with conventional photographic methods we are apparently imable to reconcile 
low distortion with high efficiency storage. Faced with this impasse, let us review 
what we mean by distortion. Since lens distortions (i.e. , barrel and pincushion) are 
stationary, they could, at least in principle, be compensated by an inversely distorted 
master print. On the other hand, it is difficult to conceive how a large number of 
various enlarger-printer devices could be adequately corrected. The second type of 
distortions are those introduced by flexible base photosensitive material, and these 
are nonstationary and present an even less tractable problem. These distortions are 
random but not completely isotropic since they generally have a strong dependence on 
the transport mechamsms used durir^ pre-ejq)0sure, exposure, and processing. In 
addition to limiting straightforward mensurational accuracy, film distortions also 
upset the correlation techniques that are quite fundamental to photointerpretation. 

The approach taken in this study is to avoid optical distortions by recording and 
retrieving at IT. Film distortions cannot be avoided, but at least with this approach 
they are not amplified by a reduction-enlargement ratio. Furthermore, when spectral 
zones are eiqjosed by multiple imaging methods, random film distortions present no 
difficulty in retrieval registration because local film distortions are identical for all 
channels. 

Photo-optical multiplexing has no complete analog in the electrical signal domain, 
but it can be compared to multichannel telemetry where many subcarriers modulated 
by Gie output of various sensors are mixed and transmitted over one high frequency 
carrier. The subcarriers are separated upon detection either by means of their 
unique frequency or by their unique phase. 
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Optically we use spatial carriers of either unique frequency or phase, but unlike 
the electrical analog we can vary the carrier frequency and phase by a simple rotation 
or translation in the plane of the carrier. It may also be helpful to think of a photo- 
graphic transparency as performing some of the functions of the high frequency carrier 
of multichaimel telemetry. Just as one would not use a subcarrier of frequency com- 
parable to the main carrier, so we are constrained to use spatial carriers that are 

3 4 5 

well resolved by an emulsion. The sampling theorem ’ and aliasing constraints 

must also be observed m optical multiplexing if no image information is to be lost 

or distorted in crosstalk. 

SUMMARY OF WORK 

Recording, processing, and retrieving are the three major task areas associated 
with the primary study objective. For recording we constructed a 1:1 multiple expo- 
sure copy camera from existing bench hardware and prepared two sets of twelve input 
targets. One set was intended for subjective evaluation and Initially consisted of 
astronaut and moon landing slides. Later, because of the poor photographic quality 
of these rather cheap slides, we switched to higher quality travel slides of the 
Washington, D.C. area. The second set of targets was designed specifically for 
objective measurements and consisted of high and low contrast resolution targets, 
step tablets, a mensuration grid, and general scenes. The dynamic range of all the 
targets was controlled so that it did not exceed a 32:1 brightness ratio, and had an 
average brightness ratio of about 20:1. A special target was also designed for cross- 
talk measurements. Nine different multiple image recordings were produced; seven- 
nine-, and twelve-channel e^qiosures of the subjective, objective, and crosstalk targets. 

Photographic processing depended on film selection. In the course of this study 
three different films were evaluated Eastman Kodak SO-343, 'Minicard, and Recordak 
AHU Type 5460. Direct reversal processing for a gamma of -2 was worked out for 
each film. In addition, a print reversal process was devised for SO-343. Since 
direct reversal of Type 5460 microfilm proved to be optimum, it was used for all the 
final multiple image recordings. 
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Liitially, demodulation and image retrieval were performed on a general purpose 
spatial filtering optical bench. Later a special 1:1 projection system was designed 
and set up, and it was used to record or measure all retrieval imagery and documen- 
tation results. 

The secondary objective of this study required the generation of synthetic target 
signals, printout of the targets, and retrieval of the targets in color. Before the first 
two tasks could be completed, however, the photofacsimile device had to be reacti- 
vated and calibrated. Although these requirements were beyond the scope of the study, 
they influenced the schedule and final results . 

The first step in generating the synthetic target signals was to specify graphically 
the carrier-modulated density distributions for a six-color bar chart, a seven-step 
gray scale, and a six-element resolution target. The density distributions were then 
translated to command language for the photofacsimile device. In turn the commands 
were reduced by a computer routine to a pimched tape. The pimched tape was, in 
effect, our simulation of the three-zone imagery stored with temporal registration in 
a buffer. 

The tapes were then used to control printout of the three targets in the photofacsimile 
device. The resulting black-and-white transparencies were then reconstructed in color 
on an optical bench decoder system. The six-bar color chart (red, green, blue, cyan, 
magenta, yellow) was photographed with Kodacolor and printed. The gray scale and 
resolution targets, although made up of additive primary components, were retrieved 
as black-and-white. The resolution of this target was limited to 4 line-pair s/mm 
by the stepping rate of the photo facsimile device. 

The final task of this study was to specify the scanner-modulator requirements for 
encoding and printing three-channel imagery at high speeds for full-format 
10 line-pairs/mm inputs. 

SUMMARY OF RESULTS 

Initial retrieval results of six- and seven-channel recordings were disappointing 
because of excessive crosstalk, noise, and limited dynamic range, hi retrospect, it 
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is clear that these problems were due chiefly to our first choice of films, SO-343. 

We were motivated to use this film because its thick, Estar base held a potential for 
outstanding dimensional stability, and its high resolution emulsion apparently offered 
sufficient storage capacity. Ejsperimentally, however, the Estar base proved to be a 
source of optical noise, and its storage capacity for carrier modulated Imagery was 
inferior to that of Type 5460 microfilm. 

Seven-channel recordings of the subjective targets on Recordak AHU Type 5460 
microfilm exhibited virtually no crosstalk and had very low noise backgrounds. Nine- 
and twelve-channel recordings had barely detectable crosstalk in two to three retrieved 
channels. Objective target measurements indicated that the effective storage and re- 
trieval gamma was approximately equal to -2 for all recordings (seven, nine, and 
twelve channels ) and was linear over an input range of 10:1 in brightness, sdelding a 
100*1 output brightness range. Resolution was maintained in excess of 
10 line-pairs/mm, and dimensional distortion was less than the size of the system 
impulse response (i.e. , less than 0. 05 mm). Control target recordings indicated 
that the potential for crosstalk was greatest in those regions of the recording where 
the composite exposure level fell considerably below the average. 

Results of the secondary study objective were very much as expected despite early 
setbacks and delays in troubleshootii^ and calibrating the photofacsimile device. The 
six-bar color chart played back very saturated additive and subtractive primary colors; 
the step tablet played back neutral with the proper brightness range; and the resolution 
target played back neutral with an upper frequency of 4 line-pair s/mm, as determined 
by the 8 line-pairs/mm printing frequency of the photofacsimile device. 
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SECTION 2 


THEORETICAL DISCUSSION AND MODELS 

HISTORY 

In a broad sense, photo-optical multiplexing began in 1899 with R. W. Wood, when 
he devised an additive three-color process in which red, green, and blue separation 
positives of a color scene were copied sequentially through three gratings of different 
spacing (i.e. , spatial frequency) onto the same plate. The gratings served as en- 
coders, retaining the essential identity of the separation positives when mixed on a 
single plate. In addition, the diffractive properties of the grating modulated images 
were used to filter back the proper colors in a simple retrieval device. By 1906 the 
father-son team of Frederic and Herbert E. Ives ’ had developed a diffraction color 
process along the lines of Wood's ideas. Gratings of either different spatial frequency 
or azimuthal orientations were employed to modulate each of three color separation 
positives. The interlace scheme of Ives is illustrated in Figure 1. Bichromated 
gelatin was used to prepare these recordings by a delicate laboratory procedure. The, 
recordings, which were in essence phase transparencies, were viewed by the diffrac- 
tion chromoscope illustrated in Figure 2. A suitable limiting aperture was placed in 



Angular Code 

CnEEN 



6909033 


Figure 1. Ives' Diffraction Process of Color 
Photography (1906) 
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the chromoscope so that the diffraction 
pattern of the red image grating was 
blocked except for its red light compo- 
nent. Similarly, only the green and 
blue light components of the gratings, 
modulating the green and blue images 
respectively, were passed. 

Between 1929 and 1939, Carlo 
Bocca, probably without knowledge 
of the prior art, developed a diffraction 
process for projecting color movies. 

The process never became commer- 
cially successful, however, perhaps 
partly because of some technical difficulties and partly as a result of the success of 

triple emulsion films. The next emergence of a diffraction process color system 

11 12 13 
appears to be W. E. Glenn ’ who applied the principles to the Eidophor system 

and then to thermoplastic recording in 1959. 

The first suggestion of applying optical multiplexing to independent input targets as 

opposed to the more restricted application to the three separations of a color image 

14 15 

was made in 1962 by A. Lohmann. ‘ Implementing the methods he suggested was 
relatively impractical, but apparently for the first time the suggestion emerged that 
a number of independent images could be stored and retrieved from a single film frame. 
In 1966 R. L. Lamberts and G. C. Higgins^® demonstrated a more practical photographic 
method for recording many independent carrier -modulated exposures . However, the 
exposures with which they dealt were binary in brightness levels and represented 

17 

digital data recording as opposed to continuous-tone imagery. In 1968 F. Bestenreiner 
demonstrated the results of a carrier-modulated multiple image process that utilized 
phase carriers. As with the methods of Lamberts and Higgins, this process was also 
nonlinear and unsuited for continuous -tone imagery, and it exhibited considerable cross- 
talk. In 1968 the author introduced the idea of linear multiple image storage and for 
the first time demonstrated a continuous -tone storage and retrieval process essentially 
free of crosstalk. It is this process that forms the basis of the current study. 



Figure 2. Ives’ Diffraction Chromoscope 


8 


The evolutionary development of photo-optical multiplexing over the past seventy 

years shows no signs of diminishing, and in fact will probably continue to grow as it 

19 

receives attention from an increasing number of scientists . 

BASIC PRINCIPLES 

Spatial carrier encoding, linear photographic processing, and coherent spatial 
filtering are the basic principles that comprise the multiplexing methods of this study. 

Spatial encoding requires that the Intensity distribution of an image be periodically 
modulated prior to exposure. We do this by placing a transmission ruling in contact 
with the film emulsion and exposing the image through it. The ruling is rotated to a 
different azimuthal orientation after each ejqposure so that a unique relationship is 
established between each image and its carrier. The ruling also samples the image 
and determines the relative position of each image spectrum in the demodulation 
process. This implies that the spatial frequency and orientation of the ruling is deter- 
mined by the number of images to be recorded. Spatial encoding can be accomplished 
with incoherent illumination in a conventional camera that has been equipped with a 
rotatable ruling in its focal plane. 

Linear photographic processing is vital not only because it is the means by which 
the unique carrier-image relationship is maintained during storage, but also because 
it provides the transition mechanism between an incoherent cop5ung procedure and a 
coherent retrieval procedure. 

^atial filtering to demodulate and retrieve individual images is accomplished by 
illuminating the multiplex recording with a plane wave, erecting its Fourier spectrum 
in the filter plane, passing one of the separated spectra through the filter aperture, 
and retransforming to focus the image. Probably most insight to the multiplexing 
problem and certainly most of the carrier specifications can be realized by a carefiil 
description and analysis of the recording’s Fourier spectrum. Thus most of our 
mathematical treatment dwells on this area. 
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PHOTOGRAPHIC STORAGE 


Since the manner in which we treat the role of photographic storage is perhaps the 
single most distinctive feature of the current multiplexing method as compared with 
other proposed schemes, reviewing the fundamentals of film storage will be worthwhile. 

Input exposure is conventionally related to output density of a processed film by a 
characteristic curve referred to as the H and D (Hurter and Driffield) or the D log E 
(Density versus log^o ^^osure) curve. An idealized D log E curve is shown in 
Figure 3 . For exposures restricted to the linear region, density is related to log 
exposure by the slope or gamma of the curve as 

D = -y log E . . (1) 



7106038 

Figure 3. Idealized D log E Curve 

I 

Of more interest to us is the relation between exposure and output transmittance , 
Density is defined in terms of intensity transmittance?^ by D = -log (5^), and the 
modulus of the amplitude transmittance T is simply the square root of ^ Assuming 
that any phase components in the emulsion are rendered negligible, amplitude trans- 
mittance of a processed film is given by 

T = E*">'/2 . (O' 
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Thus by processing the film so that 7 = -2 we obtain the required linear relation 
between input exposure (which is accomplished under normal incoherent copsring con- 
ditions) and ou^ut amplitude transmittance (which assumes an essentially coherent 
plane wave illumination) . 

Of course exposures must always be restricted to the interval between log 
and log E „ which constitutes the available input dynamic range, DR^ =log (E /E . ). 

JIXqa J.11 XiIcLa XilXiX 

If the film is to be exposed M times to different images but of the same dynamic range, 
log (E' _„/E' . ), where the primes indicate image exposures, then we have the re- 

mill 

striction given in Eq. (3) on the brightness range of each exposure: 



For example if a particular 7 = -2 process has a = 15 dB and M = 6, then 
log (E '^^^ /E ^.^ ) < 7 . 8 dB or the brightness range of each exposure should not exceed 
6:1, and for M = 12 the ratio drops to 2. 6:1, Although at first sight the restriction 
implied in Eq. (3) appears to be quite severe, two experimental conditions tend to 
soften its impact. First, for a series of random scene exposures, all maxima and 
minima eiqiosure regions are not likely to be coincident and in that respect Eq. (3) 
represents a worst possible situation. Second, the 7 = -2 requirement for linear 
storage (where coherent retrieval is anticipated) implies that the output brightness 
distribution will have a quadratic distortion. To compensate for this systematic dis- 
tortion, the original image should be compressed by an intermediate storage at 
7' = -1/2. If this is done, then the actual input brightness ratios for the six and twelve 
exposure examples become 36:1 and 6.8:1 respectively. The process we are suggest- 
ing is summarized by the following equations . 


-|-7' 


E(x,y) =S^(x,y) 

(4) 

r -i-y/2 


g7^(x,y) = T(x,y) 

(5) 
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r yy'^-y/^ 


TDmod<^-y> 

o 

{ 

E(x,y) 1 \ 

J 


(6) 


where E(x,y) is one of the original image exposures processed with y' to an inter- 
mediate transparency, ^(x,y), which in turn is exposed and processed with y to a 
recording having amplitude transmittance T(x, y). When demodulated, the recording 
yields an intensity distribution y)» which since y is fixed at -2 will be linearly 

proportional to E(x,y) only when y' equals -1/2. 

Up to this point we have only considered direct reversal processing since that was 
the type used most successfiilly in this study. In principle, print reversal processing 
could also be used to achieve linear storage. Print reversal is represented in Eq. (7). 


T = 






(7) 


where y‘ characterizes the internegative process andy is the slope of the second or 
release positive material. The storage requirement is that y 'y equal 2 and of course 
that eiqiosures are restricted to linear regions of the respective D log E curves. 

We have kept our treatment of photographic storage elementary and yet have included 
all the parameters necessary for a basic understanding of its role in this multiplex 
method. A number of more subtle effects that indicate that gamma is not always the 
constant we would like it to be can be found in the literature. Higgins has recently 
shown how the Eberhart effect can render y a function of spatial frequency; the Clayden 
and intermittency effects might tend to make y a function of exposure level when se- 
quential eiqjosures are involved; light Intensification could produce nonuniform emul- 
sion sensitivity. Knowing that these potential hazards exist is useful, but they certainly 
do not prevent us from postulating a simple imified process model. 

GENERALIZED OVERLAY MODEL 

The most direct method of introducing our functional notation is to describe the 

23 

retrieval projector as sketched in Figure 4. Optical Fourier transformation is the 
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DEMODULATED IMAGE 


T(x',y'l@F(x‘,y’) 



< 7106039 

Figure 4, Coherent Spatial Filtering and Retrieval System 


fundamental operation performed repeatedly in this projector. Lens produces the 
Fourier transform of the amplitude distribution of the (fx' ^' ) plane in the (x, y) plane; 
lens Lg transforms the (x,y) functions into the (jx, plane; lens likewise trans- 
forms the (pi, distribution into the (x' ,/ ) plane. A quasi -monochromatic point 
source of wavelength X represented by the delta function 5(pi' , I') produces a plane 
wave in the (x,y) plane as indicated by Eq. (8). 
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Of course the plane wave is truncated by the size of the collimator lens The 
projector film gate is located in the (x,y) plane; therefore when a transparency is 
inserted, its spectrum is displayed by the field lens in the filter plane as shown 
in Eq. (9). 


T(m, I) 




T(x,y) e 


27ri 

f2X 


(MX+ 4y) 


dx dy 


( 9 ) 


A spatial filter represented by the transmittance function F(^, is introduced to 
modify the spectrum, and the product T(jit, 4) F(fi, |) is transformed by the projection 
lens Lg to form a demodulated image in the (x'jy') retrieval plane. This process is 
described in Eq. (10) where the asterisk symbolizes the convolution integral. 


4) F(fi, i) 




TQjl, 4) F(n, e 




dju d^ 




T(o-, 7 }) F (x' - ff, y' - Tj) da di] 


( 10 ) 


= T(x',y) © F(x',/) . 


The focal length of the transform lens and mean wavelength \ of the illuminant 
are constants for a given projector, and these parameters determine the spatial fre- 
quency scale in the filter plane. As an example we will use Eq. (9) to calculate the 

size of the spectrum for a transparency with maximum spatial frequency w (i.e. , a 

s 

transparency that is band-limited to 2 o) ). Since the outer boundary of the spectrum 

s 

will be determined by the highest frequencies, we need only consider a typical compo- 
nent, cos 2jt w (ckx + /3y), where ot = cos 0, /3 = sin 0, and 0 designates the orientation 
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of the frequency component (e.g. , 0 = 0 for a cosine fringe varying only in the 
X direction). Rewriting the cosine as 1/2 [ e 27 riWs(®x+ py) ^ ^ 27riWg(o!x+ /3y)j 

takh^ the transform, we obtain 



_oo — oo 



The delta functions of Eq. (11) indicate that the greatest radial dimension p_ of the 

s 

spectrum is reached at the coordinate positions jj. = oico f \ and ^ = jScu f„\, but p 

9 2 1/2 ^ ^ 

equals so that we obtain the result given in Eq. (12). 

Ps' . (12) 

Thus allowing for random orientations of the w components, we know that a circular 

region of radius p would be occupied by the frequency spectrum, 
s 

Let us now examine what happens if we carrier-modulate this band-limited trans- 
parency. As stated previously, our spatial carriers are rulings and any one -dimen- 
sional ruling can be described in terms of a Fourier series as in Eq. (13). 

27 rin&) (orx+ jSy) 

P(x,y) = 2y ^n ® 

n=-<=o 
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where n's are integers including zero, a^*s are the Fourier coefficients, the 
fundamental ruling frequency, and a and /S have the same meaning as when previously 
used. Even though the ruling of Eq. (13) extends over a two-dimensional plane, it is 
one-dimensional in the sense that it is described by a single set of Fourier coefficients: 


2W Q! 2i0 B 
,0 o'^ 


= I I 


P(x,y) e 


-27rinco (Q!x+ /3y) 


dx dy . 


2<i) a 2o) B 
o o'^ 


(14) 


Other than fundamental frequency, the a^^'s specify all the physical properties of the 
ruling such as maximum and minimum transmittance, duty cycle, and slit shape. 

The Fourier spectrum of this ruling as it would be formed in the retrieval projector 
is given in Eq, (15). 


I) 






o’ fgX 





(15) 


By sandwiching the ruling and the band-limited transparency together in the film 
gate and describing the filter plane distribution, we can gain some insight to sampling 
requirements. This convolution process is given in Eq. (16). 

T(x,y) P(x,y) = T(m,|) ® P(p, ^) 

00 

= ^ a^T(M,4) ® -nofW^, ^ -n/?co^^ . (16) 

n=-oo 

From our previous calcvilations we know that the diameter of T(fz, ^) is 2w f„X and 

s ^ 

that the distance between centers for adjacent delta functions is w f_X. Figure 5 sketches 

o ^ 
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some of the spectrum characteristics indicated in Eq. (16) for the cases of 

2w > Cl) , 2co = 0 ) and 2Ci) < ct) , which from the graphical representation are quite 
SOSO so 

clearly the cases of under sampling, just adequate, and oversampling respectively. 
Overlap of the replicated transparency spectra is symptomatic of aliasing in an vaider- 
sampled signal. Oversampling, on the other hand, simply spreads out the spectra in 
the filter plane and produces a greater guardband between them. 


211)5 
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Figure 6 . Geometry for an M Channel 
Spectrum 


From the Figure 5 sketch we are in a 
good position to specify the minimum 
carrier frequency and angular rotation 
increments required to accommodate an 
M channel recording (where M is num- 
ber of channels) in terms of the bandlimit 
of the images such that there will be no 
overlap of the first order spectra. Since 
the amplitude carriers produced a sym- 
metric replication of spectra, only half 
the frequency plane is available for cod- 
ing. The situation is sketched in Figure 
6 where the length of the total available 
bandwidth as indicated by, the shaded 
semi-annulus is to)^. Assigning equal 
and contiguous segments of this semi- 
annulus to each channel leads to the 
requirement that 


A0 



(17) 


When M is greater than 4 or 5, then to a good approximation we can assume that 
A 00 ) equals 2co and so the carrier frequency requirement is obtained as shown in 
Eq. (18). 


0 ) = — 20 ) . (18) 
O IT S 


18 



To describe the carrier rotations, we simply introduce the summing index 
m = 1, 2, M into Eq. (13) to obtain Eq. (19): 


P (x.y) 
m' 


- 1 


27rinctj 






n = -“ 


( 19 ) 


where 


oi = cos 
m 


9+ (m - 1) 


7T 

M 


and 


/3. 


m 



e + 



The generalized overlay M channel recording is then described by Eq. (20), 


T(x,y) 



T-y/2 




( 20 ) 


^ th. 

where E^(x, y) is the latent exposure of the m image. When linear photographic 
storage is achieved, the brackets of Eq. (20) are removed along with at least the 
mathematical possibility of generating crossproducts. In fact, any mechanism in the 
storage process that would make it unrealistic to describe the recording as in Eq. (21), 


T(x,y) = 


■ M 
m=l 


E (x, y) 
m' 




( 21 ) 


where and Kg could be additive and multiplicative constants, would imply that the 
imique relationship between image and carrier had been lost and some degree of 
Crosstalk should be expected. 

* 

Considering Ej^(x, y) as the latent exposure obviates the questions of spectral 
sensitivity and exposure duration since this quantity is the result of integrating over 
time and wavelen^h. 
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th 


To demodulate the m channel we first Fourier transform Eq . (21) (let Kj = 0 
and Kg = 1) as shown in Eq. (22). 


M 


T(fi, ^) = ) } a E ( 4^ - noi CO , -p^r - nB co 

Z_y /-/ n m y f 2 >. m o’ f 2 \ ^m o 


m=l n=-«> 


( 22 ) 


We must then introduce a spatial filter to pass one of the m^^ chaimel’s fundamental 
orders. If we define the transmittance of a clear circular aperture of diameter 2p 

s 

in an opaque screen as 

1 + r) < p„ 


(M.ll 2p^) 


1/2 


0 VM + ? >/ > 


(23) 


th . 


then the retrieval filter for the m image order) is given by Eq. '^24). 


9 = (t* - ' ^m“o¥ ) 


(24) 


The amplitude distribution transmitted throi^h the filter is 


T(K. 5) F„(C, 0 = - VoV > « - '’m“oV) 


(26) 


which is then transformed by lens L„ of the projector to yield 


T(x',y) ® F^(x',y) = a^E^(x',y^) e 


27 rii‘^ f- / CK^x' + fi y' \ 
o f« m ) 


(26) 
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However, only the modulus squared of Eq. (?.6) is detected as the demodulated image 
shown in Eq. (27). 

1 T(x',y^) ® F^(x',y')|2 = I aj^ E^(x,y) . (27) 

The brightness distribution of the retrieved image is quadratic ally distorted with 
respect to the input exposure as predicted during our discussion of photographic 
storage. This distortion is easily compensated, however, by putting the actual image 
through an intermediate "square root" process before exposure in the multiplexed 
recording. The overall brightness level of the retrieved image is weighted by | a^^ [ 
methods of optimizing these fundamental Fourier coefficients are treated in the 
next topic. 

In addition to establishing the basic feasibility of a process, an important function 
of any theoretical model is to aid in specifying the conditions imder which the process 
should be conducted. If the model is too general, then specific physical parameters 
may be omitted; if the model is too specific, it may become unwieldy and still not 
furnish enough information about the general process. Our hope is that this treatment 
has struck a balance between these extremes. 

ENCODER OPTIMIZATION 

3h the previous section we assumed that demodulation would always be accomplished 
with one of the fundamental orders and thus retrieval brightness would be proportional 
to I I To obtain the most favorable signal-to-noise ratio, we wish to maximize 
I a^^ 1 2, which is a measure of modulation efficiency. 

The simplest carrier that we can specify is the sinosuid of a single harmonic fre- 
quency (i.e. , a periodic function where all the a^'s are 0 for n greater than 1) as 
given in Eq. (28). 


P(x, y) = Uq + a^ 


2irico (o!x+ jSy) 
e -I- a 


-2iriw^(ax+ j3y) 
e 


(28) 
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If we choose the origin of P(x,y) so that it is an even function, then a^ = a_^, and 
Eq, (28) is simply a cosine fringe as shown in Eq, (29). 

P(x,y) = + 2a^ cos 2irw^ (Q!x+ j5y) . (29) 

The problem is to specify a^ and to maximize fringe visibility. The first reaction 
is to let a^ go to zero, but this would only lead to a frequency doubling since film is 
capable of storing only the modulus and not the phase information of the carriers. 
Thus we have the condition that a^ > and so the highest contrast sinusoid 
carrier is 


P(x, y) 


1 

2 


1 + cos (ax + /3y) 


(30) 


where a^ = 1/2 and a^^ = 1/4. 

If a transparency with uniform modulation as specified in Eq. (30) is placed in a 
retrieval projector, 25 percent of the film gate energy will go to the zeroth order, 

6. 25 percent will be diffracted into each of the fundamental orders, and the remaining 
62. 5 percent will be absorbed by the transparency. 

Now let us consider the carrier having the largest number of harmonic components, 
namely a rectangular ruling as shown in Figure 7, where P^ and P^ are the maximum 


P (x.y) 



Figure 7 . Rectangular Ruling Transmittances 
(a = 1, ^ = 0) 
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and minimiun transmittance levels respectively, d is the slit width, and dw^ is the 
duly cycle. Using the half range of Eq. (14), since the function is even, we calculate 
the Fourier coefficients as shown in Eq. (31). 




d/2 




cos 27 mcii X dx + P 
o 


l/2co^ 

» I 

d/2 


cos 27mw xdx 
o 


■ = Pq sine n + ^ - P^ j doj^ sine ndco^ , 


(31) 


where sine z = (sin ttz/ttz). 

From Eq. (31) we see that a^ = and = (P^ - P^) dtu^ sine dw^. First, it 

is clear that P^ should go to zero to maximize a^ for any duty cycle. Second, by 
rewriting as 

_ Trdw^ , (32) 

1 IT o 

it can be seen by inspection that a^ will be maximum for a 50 percent duty cycle (i.e. , 
dcj^ = 1/2). Normalizing the peak transmittance level (P^ = 1), making the ruling bars 
opaque (Pq == 0), and using a 50 percent duty cycle 3 i,elds a maximum modulation 
efficiency of | a^ | ^ = (1/ir)^ *=“10 percent for a square wave ruling, which makes it 
about 60 percent more efficient than the highest contrast sinusoid. However, just as 
with the sinusoid, 25 percent of the film gate energy goes to the zeroth order, but 
only 50 percent is absorbed by the transparency and the remaining 5 percent is dis- 
tributed among all the higher order harmonics. 

Since within the class of real, non-negative going functions, the only other possible 
carriers fall somewhere between the two extreme examples we have treated, we con- 
clude that a high contrast square wave ruling provides the optimum encoder. 
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SQUARE ARRAY INTERLACE MODEL 

We develop here the theory behind the secondary study objective and continue as 
much of the previous functional notation as is practical. Based on the principle of 
optical subtraction, we had proposed that buffered three-zone spectral imagery 
could be printed out in a square array geometry so that registration of the three zones 
was preserved and color information could be retrieved. The printout geometry is 
shown in Figure 8, where each sample element is composed of four subsamples 
(A, B, C, D). By defining a high-contrast square wave ruling in the x direction as 

CO 

27rinco X 

P(x) = 2^ a^ e , (33) 

n = -oo 

with 


y 





7 I 06043 

Figure 8. Square Array Geometry 
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and in the y direction as 


P(y) = 


2irU(oy 

o*^ 


_g = _oo 


1 . i, 

2 sine 2 , 


we can describe the square array transmittance by superposition as 


T(x,y) = AP(x) P(y) + BP(x) P(^y + ^ ^ 


+ CPf x + 


^)pW.Dp(x.^)p(y.^ 


and the Fourier spectrum is 


OO 00 


T(4, ^ 2 V£ « (f* - « - ^“o'2^ 


n=_eo ^ = _C 


A+ B C e“’^'^+ D e 
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At various orders as designated by the n, & indices, we obtain the following sums 
and differences: 


n = o = a^6(f^,|)[A+ B + C+ D] (39) 

i = 0 

i = 0 

TO^»^)|n=0 = Vl ® 

& = ±l 

T(lJi, ^)L =±1 = aj a (m + ^"^"o^ 2 ^) [ A - B - C + D] ( 42 ) 


If we make D = A + B + C, then the subsample spectra separate in a very orderly 
maimer as indicated in Eqs. (43) through (46). 


T(m, 0 = 2aQ6(lLt,^)D (43) 

& = 0 

T(ju,|) = -2a^aQ M B (44) 

i = o 

T(m, 4) j^=o " 5 (^li, I +^0^2^) C (45) 

^=±1 

^>|n=±l " ^ ^2^’ ^ ^2^) ^ 

^=±1 
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This distribution is shown in Figure 9. An identical distribution can be obtained by 
letting A = B' + C', B = C' + A', C = B' + A', and D = O in Eqs. (39) through (42). The 
first scheme we refer to as additive because for normal color imagery A, B, and C 
would be the additive primary separations (red, green, and blue) and D would be the 
integrated full-color image as recorded on a panchromatic black-and-white film. The 
second scheme is subtractive since the subsample exposures would be equivalent to 
cyan (white minus red), magenta (white minus green), and yellow (white minus blue). 
However, whether the three zones to be recorded are in the visible, ultraviolet, or 
infrared does not matter. As long as the logic of one of these two schemes is followed, 
we can store the information in registration and subsequently obtain spatially sepa- 
rated zone spectra. By inserting color filters over the various zone spectra ^ all 
channels can be projected simultaneously in perfect registration to produce true or 
false color images. 
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Figure 9. Square Array Fourier Spec- 
trum Geometry 
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SECTION 3 


EXPERIMENTAL APPARATUS, PROCEDURES, AND RESULTS 

Only the overlay multiplexing experiments are described in this section; results of 
the square array interlace experiments involving the photofacsimile device are given 
in Section 4. 

GENERAL 

To record twelve channels of 10 line-pairs/mm imagery, Eq. (18) indicates that 
the encoding ruling must be at least 80 ^p/mm. - Riitially we used a carrier of that 
frequency, but all the final recordings reported in this study were prepared with a 
100 ^p/mm Ronchi ruling, which provided a small guardband even for the twelve- 
channel experiments. Encoding angular mcrements of 15, 20 and 24.8 degrees were 
.used for the twelve-, nine- and seven-channel exposures respectively. Six-channel 
experiments were discontinued early in the study because they did not permit storage 
of all seven ERTS separations on a single film frame. Twelve- and nine-channel re- 
cordings provided a limited means for documenting system parameters as a function 
of the number of images multiplexed. 

COPY CAMERA 

Photographs of the copy camera used for all of the overlay recordings are shown in 
Figure 10. The side view (Figure 10a) shows the major camera components: an 
Aristo Trans-Luminator, Model T-16, Series 3200 served as an extended uniform 
incoherent source; a Leitz 2 x 2 in. slide transport was used for holding the 35 mm 
object transparencies; a Fax-Nikkor f/5. 6 (160 mm) was used for the copy lens; and 
a simple ruling holder-pressure pad fixture provided the film plane. The rear view 
(Figure 10b) shows this fixture in more detail. The Ronchi ruling is a 2 x 2 in. High 
Resolution Plate copy of a 100 ^p/mm metal master, and it is cemented to a brass 
ring that can be rotated in the square aluminum plate. Studs on the back side of the 
ring enable rotations to be indexed in total darkness. Pins on the front side of fixture 
plate are used to locate and hold 35 mm film chips by means of their perforation holes. 
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C 1:1 COPY LENS 
D FILM PLANE 


A LOCATING PINS 

B ROTATABLE 
RONCHI RULING 

C PRESSURE PAD 




b. Rear View 
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Figure 10. Copy Camera 
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A pressure pad is mounted on a second aluminum plate that is hinged to the first on 
the right side. During exposure the second plate is closed and locked by means of 
three thumbscrews, and the spring-loaded pressure pad is actuated. 

The distance between the object and image plane is 639.5 mm, and the lens is 
focused at a location just midway between the two planes. No special alignment pro- 
cedure was used in setting up this camera. However, since all the machined surfaces 
of the optical bench carriers and uprights are fairly well squared off, they served as 
good, but by no means precise, references. For extremely precise alignments the 
three-point mounting feature of these uprights can be utilized to set two planes parallel 
by autocollimation, and the lens axis can be made perpendicular to the planes by the 
Boys point method. It was not necessary to use such techniques for this study. 

TARGET PREPARATION 

Two sets of targets were prepared for inputs: subjective targets, which in the final 
recordings consisted of the 35 mm travel slides shown in Figure 11, and objective 
targets as shown in Figure 12. The subjective targets were black-and-white copies 
made from color slides with reversal processed Panatomic-X film that was exposed 
so the maximum density difference would not exceed 1.5 and the average density dif- 
ference would be about 1.3. These transparencies therefore had a considerably 
greater brightness range than was indicated desirable by the theoretical model in 
Eq. (3). However, at the time the targets had to be prepared we decided it was better 
to risk the possibility of some crosstalk rather than overly restrict the input image 
contrast. In retrospect it is clear that an extended range reversal process equivalent 
to 7 = -1/2 storage would have been preferable, but additional photographic investiga- 
tions would have been required to work out such a process initially. 

Only five objective targets were prepared as such, and the additional seven re- 
quired to make up to twelve exposures were selected from the subjective set. The 
three resolution targets had the standard spatial frequency progression given in 
Table 1. The high-contrast positive target had a density difference of AD = 1.3, the 
high-contrast negative a AD = 1.10, and the low-contrast positive a AD = 0.4. 
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Figure 11. Subjective Input Targets 

Seven-Channel: a through g 

Nine-Channel: a through i 

Twelve-Channel: a through 1 
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Figure 12. Objective Input Targets 

Seven-Channel: these targets plus a and g of Figure 11 
Nine-Channel: these targets plus a, g, h, j of Figure 11 
Twelve-Channel: these targets plus a, g, h, j, 1, e of Figure 11 


The step tablet density calibration as measured on a Macbeth Densitometer, Model 
TD-102 is given in Table 2. This target was made by overlaying the first seven steps 
of a Kodak step tablet as given in the first row with a strip of uniformly fogged film of 
0. 08 density to form the second row and a neutral density filter strip of 1. 0 to form 
the third row. 

Calibration of the mensuration grid target is given in Table 3 and the line-colunm 
designations are shown in Figure 13. This target was calibrated by averaging three 
measurements for each line and column made on a Scientific Instruments comparator 
that had a smallest division scale of 0. 01 mm. 

The final target designed specifically for crosstalk measurements was simply a 
0. 9 neutral density filter placed over one half the 35 mm field and rotated with each 
exposure about the optical axis to produce the fan-shaped ejqjosure shown in Figure 14. 
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Table 1. Resolving Power Target ^ Progression 
(Values for Individual Targets) 


Group Element RP 

No. No. £p/mm 


Group Element RP 

No. No. /!p/mm 

1 0.250 

2 .280 

2 3 .315 

4 .353 

5 .397 

6 .445 


1 16.0 

2 17.95 

. 3 20.16 

4 22.62 

5 25.39 

6 28.51 

1 0.500 

2 .561 

, 3 .630 

4 .707 

5 .793 

6 .891 


1 32.0 

2 36.0 

c 3 40.3 

4 45.3 

5 50.8 

6 57.0 

1 1.000 

2 1.12 

0 3 1.26 

4 1.41 

5 1.59 

6 1.78 


1 64.0 

2 71.8 

„ 3 80.6 

^ 4 90.5 

5 102. 

6 114. 

1 2.00 

2 2.24 

, 3 2.52 

4 2.83 

5 3.17 

6 3.56 


1 128. 

2 144. 

7 3 161. 

4 181. 

5 203. 

6 228. 

1 4.00 

2 4.49 

„ 3 5.04 

4 5.66 

5 6.35 

6 7.13 


1 256. 

2 287.4 

„ 3 322.6 

® 4 362.0 

5 406.4 

6 456.2 

1 8.00 

2 8.98 

„ 3 10.1 

4 11.3 

5 12.7 

6 14.3 
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Table 2. Step Tablet Calibration 


Step 

Density Measurement 

Row 1 

Row 2 

Row 3 

1 

0.06 

0.14 

1.06 

2 

.20 

.28 

1.20 

3 

.34 

.42 

1.34 

4 

.48 

.56 

1.48 

5 

.62 

.70 

1.62 

6 

.76 

.84 

1.76 

7 

.90 

.98 

1.90 


Table 3. Mensuration Grid Calibration 


Measurement Average (mm) 

Lines 

Columns 

Lj = 31.26 

C, = 15.62 C. = 15.63 

1 5 

Lg = 31.27 

Cg = 15.62 Cg = 15.63 

Lg = 31.27 

Cg = 15.63 C.^ = 15.63 

L. = 31.26 
4 

C^ = 15.63 C„ = 15.64 

4 o 

L, = 31.26 


5 




•-4 


JU :: !_L5^! !! 2! JL 

7 1 06046 

Figure 13. Line -Column Designations for 
Mensuration Grid 
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Figure 14. Crosstalk Control Exposure Pattern 

a. Seven -Channel 

b. Nine-Channel 

c . Twelve -C hannel 

MULTIPLEXED RECORDINGS 

The procedure followed in making the multiple image recordings was partially 
indicated in the copy camera description. Film handling for all the final recordings 
had to be accomplished in total darkness since AHU 5460 is panchromatic. However, 
infrared light and a Snooperscope were used to facilitate certain procedures. 

One of the first recording defects noticed was Newton rings caused by interference 
of light reflected in the air cavity between the AHU emulsion and the Ronchi ruling 
emulsion. The best way to suppress this interference phenomenon is to anti-reflection 
coat the ruling surface and thereby eliminate the resonant cavity that supports the 
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NOT REPRODUCIBLE 


rings. Even though high efficiency anti-reflection coatings are simply not practical 
for application to a gelatin surface, they could be evaporated over or under a metallic 
ruling copied onto a glass substrate. The lead time required to fabricate such a railing 
(8 to 10 weeks) precluded this approach in the present study. 

An alternate method of eliminating the reflection is to fill the cavity with a liquid 
having a refractive index matched to the gelatin. This was the approach we took, and 
it required the application of a few drops of xylene between the ruling and emulsion 
before each exposure. The Snooperscope was required to execute this operation. It 
must be emphasized that this liquid gating procedure was simply an expedient in these 
experiments and that an anti-reflection coated ruling would have to be used in any 
practical system. 

All exposures were made at a full lens aperture of f/5. 6 at an average exposure 
time of about 1/5 of a second each, although some targets required as little as 1/10 
of second and others as much as 1/2 of a second for AHU 5460 film. Photographs of 
the multiplexed recordings are shown in Rgure 15, and details of the modulation pat- 
terns are evident in the photomicrographs in Figure 16. 
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Figure 15. Multiplexed Recordings 

a. Seven-Channel Subjective d. Seven-Channel Objective 

b. Nine-Channel Subjective e. Nine-Channel Objective 

c. Twelve -Channel Subjective f. Twelve-Channel Objective 
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a. Seven-Channel 



b. Nine-Channel 



c. Twelve -Channel 

7106047 

Figure 16. Photomicrographs of Typical Seven-, 
Nine-, and Twelve-Channel 
Recordings 
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PROCESSING JUDIES 


The first photographic storage was worked out for Eastman Kodak's SO-343. It 
was a print reversal method in which the negative was processed as indicated in 
Table 4 to produce the D log E characteristics shown in Figure 17, and the positive 


Table 4. Gamma = 1 Process for SO-343 
(All solutions at 70°F) 



Step 

Duration 

1. 

Develop (Tech/Ops Extended Range 
Developer*) 

20 min 

2. 

Wash 

30 sec 

3. 

Fix (E.K. Rapid Fixer) 

2 min 

4. 

Wash 

30 sec 

5. 

Hypo, neutralize (Heico Inc. Perma Wash) 

2 min 

6. 

Wash 

2 min 

7. 

Photo Flo (E.K. 200 solution) 

30 sec 

8. 

Air dry 

15-20 min 


* 

Tech/Ops Extended Range Developer formulation 


Ingredient Amoimt 


Distilled water 

750 ml 

Potassium hydroxide 

15 g 

Potassium sulfide 

100 g 

Graphidone (Hunt Chemical Co.) 

24 g 

E.K. Antifog #1 

6g 

Potassium bromide 

4 g 

Acetic acid to pH 10 


Distilled water to make 

1000 ml 
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print was processed as specified in Table 5 to yield the D log E curve plotted in 
Figure 18. This combination of gammas satisfied the requirement of Eq. (7) that 
y'y/2 equal 2, 


Table 5. Gamma = 2 Process for SO-343 
(All solutions at 70° F) 


Step Duration 


1. 

Develop (Acufine Inc. Acufine Film Dev.) 

5 min 

2. 

Wash 

30 sec 

3. 

Fix (E.K. Rapid Fixer) 

2 min 

4. 

Wash 

30 sec 

5. 

Hypo, neutralize (Heico Inc. Perma Wash) 

2 min 

6. 

Wash 

2 min 

7. 

Photo Flo (E.K. 200 solution) 

30 sec 

8. 

Air dry 

15-20 min 


Operationally this process was very awkward because an iterative approach was 
needed to optimize both the multiple exposures and the contact print exposure. The 
multiple ejqiosures on the negative could not be optimized independently because the 
extremely fine detail of the modulation (as shown in Figure 1 6) is in the range where 
conventional microdensitometry methods become unreliable. Ultimately, print 
reversal may be the preferred storage method especially where a large number of 
positive duplicates are to be made. No doubt this approach could have been worked 
out successfully once input target exposure limits were fixed, but we abandoned 
print reversal early in this study primarily because of its relative inflexibility and 
long turn-around cycle for an experimental program. 
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To speed up the recording and storage process we next worked out a direct 
reversal process for SO-343 (as specified in Table 6) that produced the D log E 
curve shown in Figure 19, Despite the tremendous dynamic range exhibited in this 
D log E curve, SO-343 did not appear well suited for multiple imagery based on the 
observation of demodulated images. The reasons appear to be two-fold: the thick 
Estarbase has poor optical quality that is e:diibited as random noise in a dark-field 
viewing system, and resolution seems to be rapidly lost at higher exposure levels 

Table 6. Gamma = -2 Process for SO-343 
(All solutions at 68® F) 


Step Duration 


1. 

Develop (E.K. D-19) 

3 min 

2. 

Wash 

30 sec 

3. 

Bleach (E.K. R-9) 

2 min 

4. 

Wash 

30 sec 

5. 

Clear (E.K. CB-3) 

30 sec 

6. 

Wash 

30 sec 

7. 

Expose to floodlight (500 watts at 30 in.) 

30 sec 

8. 

Develop (E.K. D-95) 

1.5 min 

9. 

Wash 

30 sec 

10. 

Fix (E.K. Rapid Fixer) 

2 min 

11. 

Wash 

30 sec 

12. 

Hypo, neutralize (Heico Inc. Perma Wash) 

2 min 

13. 

Wash 

30 sec 

14. 

Photo Flo (E.K. 200 solution) 

30 sec 

15. 

Air dry 

15-20 min 
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perhaps because of its relatively poor antihalation properties. The origin of base 
noise appears to be particles of different refractive index either embedded in the sub- 
strate or in a backing coat. Normally the phase differences produced by these par- 
ticles would not be detected in diffuse or incoherent illumination, but our retrieval 
device (since it is a dark-field system, with some properties in common with a 
schlieren system) enhances the contrast of phase objects. 

The third process evaltiated durir^ this study was one we had used for several 
years both for multiple imagery and for the preparation of sequentially exposed tri- 
color images. The recording material was Recordak AHU, Type 5460, Microfile film. 
It was direct reversal processed as specified in Table 7 and yielded the D log E char- 
acteristics plotted in Figure 20. The of this process is relatively low, which 

Table 7. Gamma = -2 Process for AHU Type 5460 Microfile 
(All solutions at 68°F) 


Step Duration 

X. Develop (E.K. DK-50, 1 part stock to 2 parts water) 5 min 

2. Wash 1 min 

3. Bleach (E.K R-9) 3 min 

4. Wash 30 sec 

5. Clearing bath (E.K. CB-3) 30 sec 

6. Expose to floodlight (500 watts at 30 in.) 30 sec 

7. Develop (same as Step 1) 5 min 

8. Wash 1 min 

9. Fix (E.K. Rapid Fixer) 2 min 

10. Wash 30 sec 

11. Hypo, neutralize (Heico, Inc. Perma Wash) 2 min 

12. Wash 2 min 

13. Photo Flo (E.K. 200 solution) 30 sec 

14. Air dry 15-20 min 
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should adversely influence the dynamic range capability in multiplexed ej^osures. On 
the contrary, however, this direct reversal process produced the best retireval, both 
in terms of minimum crosstalk and maximum retrieval brightness range. 

The final process evaluated was one worked out at the Eastman Kodak Research 
Laboratory several years ago for Minicard film. This process is somewhat unique in 
that the first step calls for slightly prefixing the film, which tends to differentially 
dissolve part of the une?q)Osed silver halide more than that which is exposed. The ef- 
fect is to lower the base density and gamma in a manner similar to a prefogging expo- 
sure. Formulation of the process is given in Table 8 and the D log E characteristics 
in Figure 21. D3mamic range and D are considerably greater than for Type 5460 

Table 8. Gamma = -2 Process for Minicard Film (All solutions at 68°F) 



Step 

Duration 

1. 

Prefix (1 part fixer* to 18 parts water) 

2 min 

2. 

Wash 

1 min 

3. 

Develop (Boiler’s Developer^) 

3.5 min 

4. 

Wash 

1 min 

5. 

Bleach (E.K. R-9) 

3 min 

6. 

Wash 

1 min 

7. 

Clearing bath (E.K. CB-3) 

1 min 

8. 

Expose to floodlight (500 watts at 30 in.) 

30 sec 

9. 

Develop (E.K. D-95) 

2 min 

10. 

Wash 

30 sec 

11. 

Hypo, neutralizer (Heico, Inc. Perma Wash) 

2 min 

12. 

Wash 

2 min 

13. 

Photo Mo (E.K. 200 solution) 

30 sec 

14. 

Air dry 

15-20 min 

*Kodak 

fixing bath F-5 . 



^Boiler’s developer formulation: 

Sodium sulfite (des.) 60 g Potassium bromide. . 10 g 

Kodak Elon developing agent . 20 g Water to make 1 liter 

Phenidone (Ilford Limited) .. lOg pH = 10. 50 
Sodium hydroxide lOg 
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Figure 21, D log E Curve for Minicard Film Gamma = -2 Process (Table 8) 
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film, and since Minicard film is a higher resolution material, we would expect it to 
be better suited for multiplexed recording. Experimentally, however. Minicard 
proved slightly inferior, both in terms of crosstalk levels and its ability to retain 
shadow or low exposure detail. 

RETRIEVAL PROJECTOR 

A photograph of the retrieval projector is shown in Figure 22, and the schematic 
in Figure 23 shows the placement of the optical elements. The lamp house contains a 
100 watt Osram mercury arc, which is spectrally filtered to the 546 nanometer green 
line and focused into the pinhole aperture (1 . 0 mm diameter) . The pinhole defines a 
source plane and removes stray backgroxmd light. A slightly different film gate 
arrangement than described previously in Figure 4 is used, but the Fourier 



A LAMP HOUSE 
B PINHOLE 
C CONDENSER 
D ROTATABLE PRISM 
E MULTIPLEX RECORD 


F SPATIAL FILTER APERTURE 
G COPY LENS 
H RECORDING CAMERA AND 
RETRIEVED IMAGE 
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Figure 22. Photograph of Retrieval Projector 
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DEMODULATED 

IMAGE 



7106054 


Figure 23. Schematic Diagram of Retrieval Projector 

transform operations are identical. First, the film gate is located in a conventional 
position just beyond the second condenser lens, rather than between condenser ele- 
ments. Second, a thin prism of diopter number 100 “ 5. 5 is inserted between 

the condenser and film gate. The prism diopter or refraction angle is made just 
equal to the fundamental diffraction angle of the carrier so that by rotating the prism 
the various first order spectra are diffracted onto the optical axis. Ideally, to mini- 
mize potential copy lens aberrations, the pinhole or Fourier transform array would 
be imaged into the copy lens iris stop. However, the required aperture stop for 
spatially filtering the 10 ^p/mm spectra of the multiplexed recordings is 2ojgL^x 
3. 0 mm (see Figure 23) and the Fax-Nikkon could not be stopped down that small. 
Therefore, we placed a spatial filter screen in front of the lens rather than disas- 
semble the lens to insert the smaller iris stop. 

SUBJECTIVE TARGET RETRIEVAL RESULTS 

Photographs of the Fourier spectrum arrays as they appeared in the spatial filter 
plane are shown in Figure 24. One side order at a time was passed through the center 
of the spatial filter aperture, and the demodulated image was exposed in the 35 mm 
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a. Seven-Channel 
^ectrum 


b. Nine-Channel 
Spectrum 


c . Twelve-Channel 
^ectrum 
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Figure 24. Typical Fourier ^ectrum Arrays 


camera back onto Panatomic-X film. The film was processed to a gamma of 1 and 
the negatives printed at 1:1; the results are shown in Figures 25, 26, and 27 for the 
seven-, nine-, and twelve-channel recordings respectively. Overall, the seven- 
channel retrievals are best, but they are not markedly different from the other two 
sets. The most noticeable difference is the appearance of low -contrast, low- 
frequency moire fringes over several of the demodulated images (e.g. , g and h in 
Figure 26 and e, i, and k in Figure 27), The origin of these fringes is a slight non- 
linearity in the photographic storage. They can be suppressed, however, by select- 
ing carrier rotation increments that are not submultiples of 60 degrees as in the 
seven-channel recording. The crossproduct positions for the nine- and seven-channel 
spectra are compared in Figure 28. If there is the slightest deviation. A, from the 
nominal 20 degree rotation increments of the nine-ch ann el recording, then there will 
be a very small displacement between the primary order and the crossproduct as 
indicated in Figure 28a. With A in radians this displacement can produce fringes 
of Ao)^ frequency if both the primary and crossproduct orders are passed through 
the spatial filter ^erture. As A gets smaller, the fringe frequency simply becomes 
lower and perhaps more disturbing, but not lower in contrast. On the other hand, 
the 25.7 degree rotation increments of the seven-chaimel recording place any potential 
crossproducts on an axis just midway between primary orders as shown in Figure 28b, 







f g 


7 ) 06035 

Figure 25 . Retrievals of Seven-Channel Subjective Target Recording 

and safely out of the 2co = 20 ^p/mm pass band of the spatial filter. Thus, cross- 

s 

products can never coherently beat and form moir6 fringes in this case. The 
crossproduct constructions were shown for only one typical combination, but the 
explanation given can be extended to all other combinations and the same general 
conclusions will follow. Rotation increments should be chosen to place the potential 
crossproduct as far from the primary orders as possible and sufficient guardband 
maintained between these orders to avoid overlap. This last condition may require 
further oversampling on the part of the carrier. 
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Figure 26. Retrievals of Nine-Channel Subjective Target Recording 
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Figure 27. Retrievals of Twelve-Channel Subjective Target Recording 
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A0 = 25.7 deg 

7106056 


Figure 28. Position of Crossproduct Orders 
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OBJECTIVE TARGET RETRIEVAL RESULTS 


The effective modulation transfer characteristics of this multiplexed system were 
calculated from the retrieval results of the USAF ^ high contrast targets. Initially 
we had planned to make scanning photometer measurements of the three bar elements 
in the aerial image plane of the retrieval projector. Since that approach proved im- 
practical and unreliable for the probes we had available, gamma = 1 recordings were 
made on Panatomic-X film and the film subsequently measured on a Joyce- Loebl 
scanning microdensitometer. These density readings were then converted back to the 
maximum and minimum (1 ^^^) intensity levels of the exposure. Modulation 

was then calculated for each frequency as V = (Imax ” ^min)/^^max ^min)* 
results of these measurements and calculations are plotted in Figures 29, 30, and 31 
for the seven-, nine-, and twelve-channel recordings respectively. Ideally, if the 
encoding ruling performed as a point sampler (or as a line sampler in two dimen- 
sions) and the retrieval projector were perfectly coherent (i. e. , if the source were 
infinitely small), the transfer function would be unity out to 10 ^p/'mm and then 
fall abruptly to zero. This would not be a particularly desirable transfer function 
since it cotild introduce Gibb's phenomenon of edge ringing. Fortimately, the encod- 
ing ruling and film impulse response tend to produce average value sampling, and the 
finite size of the retrieval source makes the projector only partially coherent. These 
two conditions lead to an inherent falloff in the MTF curves. The major cause of 
visibility decay with higher frequencies, however, is not at all due to multiplexing 
principles, but rather to the MTP of the copy camera lens. Camera lenses can be 
designed that have reasonably flat transfer characteristics out to a particular frequency 
and then drop very rapidly. Such lenses were outside the scope of this study, but they 
are used in the television industry where large apertures are required to gather as 
much light as possible, but where high spatial frequency components could cause 
aliasing. 

There is no direct reason why the MTF should fall with an increasing number of 
channels, but it is intuitively obvious that a greater demand is placed on the film's 
storage capacity, which in turn should lead to lower signal -to-noise ratios in retrieval. 
We should also expect that with the larger number of exposures, crosstalk (if it exists) 
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Figure 29. Effective Modulation Transfer Function for 
Seven-Channel Recording 
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SPATIAL FREQUENCY IN ip/mm 
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Figure 30. Effective Modulation Transfer Function for 
Nine-Channel Recording 
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7106059 

Figure 31. Effective Modiilation Transfer Function for 
Twelve-Channel Recording 
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will become increasingly homogeneous and eventually become equivalent to a uniform 
low-level foggit^ exposure. 

If enhancement of the higher spatial irequencies of a demodulated image was 
desirable in any particular application, the required operation is quite straightforward. 
Instead of using a clear, uniformly transparent apeiiure as the spatial filter in the 
retrieval projector, a shaded or inverse filter can be placed over the clear aperture. 
This filter is designed to attenuate low frequency components most strongly, inter- 
mediate frequencies only slightly, and high frequencies not at all . An example of 
such a filter that would make the MTF of the seven-channel recording equal to unity 
out to 10 ^p/mm is shown in Figure 32. 

The effective D log E characteristics were determined by means of the low-contrast 
step tablet described in Table 2. Hetrieval measurements were made by scanning each 
step with a large aperture photometer probe. The results are shown in Figures 33, 

34, and 35 for the seven-, nine-, and twelve-channel recordings respectively. The 
gammas of all three curves are very nearly equal to -2 as predicted theoretically, 
but with the seven- and nine-channel recordings just slightly less steep than -2, and 
the twelve-channel results a little steeper. A small variation in gamma can be ex- 
pected since the film strips were hand-processed at different times with variously 
aged chemicals. The processing variations could be further reduced by automatic 
machine processing with scheduled chemistry replenishment. The important question, 
however, is not really whether the characteristic curve is precisely equal to -2, but 
whether the curve is linear. The results achieved and plotted in Figures 33, 34, and 
35 are remarkably linear, indicatii^ that no substantial distortion need be introduced 
by this multiplexii^ method. Of course, the gamma of -2 storage itself represents a 
quadratic distortion; as ejqilained previously, however, this distortion can be antici- 
pated and completely compensated by compressing the input exposures through an 
intermediate square root (i.e. , y = 1/2) storage. Exercising this precaution, we can 
store and reproduce tonal values of 100:1 dynamic range (i.e. , 20 decibels) with no 
brightness distortions. It is interesting to note that unlike conventional sensitometry, 
we cannot generate toe and shoulder regions in the characteristic curves because 
these would require exposures that themselves would violate one of the linear storage 

requirements (i.e. , exposures must be within E and E . ). 

max min 
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MODULATION 
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Figure 32. Amplitude Transmittance of An Inverse Filter That Would Nor- 
malize the Seven-Channel MTF Out to 10 Line-pairs/mm 
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AMPLITUDE TRANSMITTANCE 


RELATIVE^ LOG INTENSITY OF RETRIEVED 
STEP TABLET (OUTPUT) 



(INPUT) 
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Figure 33 . Effective D log E Curve 
for Seven-Channel Ob- 
jective Taiget Recording 
(gamma w -2) 



RELATIVE LOG EXPOSURE 
(INPUT) 

7106062 

Figure 34. Effective D log E Curve 
for Nine-Channel Ob- 
jective Target Recording 
(gamma « -2) 
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RELATIVE LOG EXPOSURE 
(INPUT) 

7106063 

Figure 35. Effective D log E Curve 
for Twelve-Channel Ob- 
jective Target Recording 
(gamma « -2) 

Dimensional distortions were documented by means of the coarse grid target 
calibrated as specified in Table 3 . Theoretically, positional accuracy is limited only 
by the sampling period in this multiplexing method, and even that factor is no limita- 
tion if the image exposure has been properly band-limited. For example, the 
100 ip/mm carrier has a period of 10 microns and so points cgnnot be located more 
accurately than ±5 microns simply because the carrier was introduced. On the other 
hand, one should not expect that points can be located with accuracy exceeding the 
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system’s spread function; since we are dealing with 10 ^p/mm imagery, the spread 
function introduces a no mina l uncertainty of ±0. 05 mm. Retrieval results of the men- 
suration grids are summarized in Table 9. This table was prepared by first measur- 
ing the raw line and column dimensions by the same comparator procedure used to 
calibrate the grid. The magnification factors shown on the top line of data were cal- 
culated as the average deviation between input and oulput grids. In other words, be- 
tween the copy lens and projector lens there was actually a very slight magnification of 


Table 9. Line and Column Dimensions (m mm) of Mensuration Grid Targets 


Line or 
Column 

Original 
Input Target 

Seven-Channel 

Retrieval 

Nine-Channel 

Retrieval 

Twelve-Channel 

Retrieval 

Magnification Correction 


1.0000 

0.9846 

0.9814 

0.9807 


31.26 

31.27 

31.25 

31.25 

^2 

31.27 

31.28 

31.27 

31.26 

^3 

31.27 

31.28 

31.29 

31.27 

^4 

31.26 

31.29 

31.30 

31.27 

^5 

31.26 

31.29 

31.28 

31.25 


15.62 

15.64 

15.60 

15.65 

^2 

15.62 

15.65 

15.62 

15.62 

^3 

15.63 

15.64 

15.63 

15.61 

^4 

15.63 

15.63 

15.60 

15.63 


15.63 

15.63 

15.62 

15.62 


15.63 

15.62 

15.62 

15.61 


15.63 

15.62 

15.64 

15.64 


15.64 

15.63 

15.65 

15.64 
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about 1:1. 012 in going from input targets to retrievals. The magnification correction 
factors were used as constant multipliers to the raw data, and the resulting products 
were tabulated as shown. Table 9 indicates that in the overall multiplexing system 
there was no pincushion, barrel, or any other dimensional distortion. The greatest 
deviation was 0. 04 mm, which is still less than the spread function uncertainty. It is 
not surprising on the basis of the optics used that distortion was immeasurable in 
these experiments. The lenses were ssrmmetric and designed for distortionless I’l 
imaging. In addition, only a small portion of the potential lens format was utilized 
in the paraxial region where distortions would be minimum even for nonsymmetrie 
lenses. What is surprising initially is that random film shrinkage was not detected. 
However, when we recall that in this process only film chips and hand processing 
were used, it becomes reasonable to eiqiect that stretching and subsequent shrinking 
would not be as severe as in a normal film transport device in which the film is pulled 
by the sprocket holes in one predominant direction. Since the variations exhibited in 
Table 9 are, for the most part, within measurement or judgment error, we conclude 
that in this study there was no detectable dimensional distortion. 

Measurements of objective crosstalk are not as straightforward to make as those 
for resolution, gray scale, or distortion. The eiq)erimental plan was based on the 
simple proposition that, in the complete absence of crosstalk, fields of uniform expo- 
sure would be retrieved as fields of uniform intensity Variations from uniform in- 
tensity are a measure of crosstalk levels. The least complicated control target we 
could devise was a clear field divided by a 0. 9 ND filter that produced an 8*1 exposure 
difference. The choice of a 0.9 ND filter was somewhat arbitrary, but it represented 
a realistic maximum A log E. By rotating the filter along with the encoding ruling, 
the jEan-shaped exposure patterns previously shown in Figure 13 were formed. 

Taking the seven-channel recording as an example, we can see that the brightest 
wedge segment running in an almost vertical direction received the maximum com- 
posite exposure, namely 7 overlapped bright-field exposures. The segments on either 
side received 6; the next two 5, then 4, 3, 2, and 1, and finally one segment on the 
bottom received no bright-field exposure. This situation is sketched in Figure 36. 
Notice that only channels 1 and 7 have wedges var 3 dng from 1 to 7 bright-field 
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I 



Figure 36. Exposure Geometry of Seven-Channel Crosstalk Target 
(numbers inside the wedges indicate the number of 
overlapped bright-field exposures) 

exposures, whereas channel 4 has the least spread in exposure extremes varying only 
between 4 and 7 bright-field exposures. Similar exposure geometries are shown in 
Figures 37 and 38 for the nine- and twelve-channel recordings respectively. 

All channels of each target were demodulated and photographed on Panatomic-X 
film (processed to a 7 = 1 ) . Densities of each wedge segment were then measured, 
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Figure 37 . Exposure Geometry of Nine-Channel Crosstalk Target 
(numbers inside the wedges indicate the number of 
overlapped bright-field ejqiosures) 

but only for the bright-field half of each retrieval target. Densities were converted 
to relative log brightness and the average log brightness calculated for all segments 
having the same number of overlapped bright-field exposures for each of the three 
crosstalk recordings. The results of these calculations are summarized in the curves 
of Figure 39. interpreti:t^ these curves requires some explanation because the phy- 
sical quantities involved are somewhat imconventional. Figure 39a indicates that 
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Figure 38. Exposure Geometry of Twelve-Channel Crosstalk Target 
(numbers inside the wedges indicate the number of 
overlapped bright-field exposures) 


channels having, between 4 and 7 overlapped exposures will play back with relatively 
imiform intensity over the bright half field. In other words (referring to Figure 36), 
channel 4 was retrieved with most JBeld uniformity and therefore least crosstalk. 
Channels 5 and 3 where next best, followed by channels 6 and 2. Channels 7 and 1, 
however, exhibited considerable falioff in brightness for those segments of the field 
that had only one bright-field exposure. Similarly, curve 39b indicates that crosstalk 
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NUMBER OF OVERLAPPED BRIGHT- FIELD EXPOSURES 

7106067 

Figure 39. Summary of Crosstalk Measurements for 
Seven (a). Nine (b), and Twelve (e) Chan- 
nel Recordings 


was minimum for those channels that had between-4 and 9 overlapped ei^josures*. 

Thut 3 channels 4, 5, and 6 exhibited the least amount of crosstalk in the nine-channel 
recording. Finally, the curve in Figure 39c is fairly flat between 3 and 10 overlapped 
ejq)Osures, which from Figure 38 indicates that the eight channels from 3 to 10 showed 
reasonably uniform half field retrievals with very little crosstalk. It must be pointed 
out that the channel number designation corresponding to relatively low and high cross- 
talk levels is peculiar only to the crosstalk control targets and in no way is to be car- 
ried over to either the subjective or objective target eJ 5 >osure series. 

Two very significant results emerge from these crosstalk data. The most glarii^ 
deviation is the falloff in retrieval brightness for those regions of the multiplexed 
recording that have received a relatively low composite exposure level; that is, in 
terms of the Figure 39 data, those regions receiving less than 3 or 4 bright-field 
exposures . Nothing in the theoretical multiplexing model suggests this phenomenon, 
but our model relied heavily on the assumption that a photographic emulsion was 
nearly an ideal storage device. Unfortunately, a number of well-known emulsion 
effects will tend to upset our linearity conditions, especially when exposures are 
made sequentially. 

One effect that may come into play is light Intensification. The sensitivity or speed 
of an emulsion can be increased if it is subjected to a short-duration high-intensity 
fogging exposure. Of course when only two or three exposures have been overlapped, 
this type of Intensification is less likely to occur than when more exposures superim- 
pose. The Clayden and intermittency effects are also likely to introduce nonlinearities 
where few exposures are overlapped. Probably certain emulsion effects that are 
difficult to control account for the crosstalk or falloff displayed on the left-hand side 
of the Figure 39 curves. The falloff evident on the right side of the Figure 39c curve 
is probably a straightforward saturation effect. The average exposure level was 
slightly high for a twelve-channel recording and thus the film’s storage capacity was 
somewhat overtaxed beyond nine exposures. 

The second result of these measurements is that crosstalk becomes less detectable 
as the ensemble exposure becomes more uniform. Thus the more channels recorded, 
the more likely is the ensemble' exposure to approach a constant over macroscopic 
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regions (i.e. , the multiplexed recording will tend to appear more like a uniformly 
fogged film as the nxunber of overlayed exposures is increased) . 

Resolution, gray scale, and dimensional distortion data can be assimilated prac- 
tically at a glance from the objective target results. A little more effort may be 
required to understand the implications of the crosstalk data since there are no stan- 
dard or traditionally familiar methods of presenting these data in the two-dimensional 
optical domain. 
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SECTION 4 


GENERATION OF THREE ZONE TARGETS (INTERLACE EXPOSURES) 
GENERAL 

The previous multiplexing method involved the use of overlapped carriers that 
resulted in a complex pattern of exposure elements extending randomly over the 
emulsion format (e.g. , as shown in the photomicrographs of Figure 16). The method 
described in this section (the theory of which was developed previously) produces a 
geometry that is very well ordered with respect to the arrangement of sample ele- 
ments. Each sample element that is part of a precise periodic array consists of four 
subsamples having either additive or substractive interlace geometry. Because this 
multiplexing method is exactly periodic in two dimensions, the square array approach 
is well suited for printout by a modulated rotating mirror scanner, CRT, or flying 
spot scanner. 

The basic feasibility of the square array method had been demonstrated previously 
by means of dichroic encoding filters consisting of red, green, blue and, white (clear) 
subsample elements. The encoding filter was held in contact with a panchromatic film 
and a color scene exposed through the filter. Subtractive logic filters consisting of 
cyan, magenta, yellow, and black (opaque) subsample elements have also been fabri- 
cated and balanced for tricolor photography. In this study, however, encoding was 
accomplished for the first time by simulating the scanner printout geometry equivalent 
to exposure through a subtractive square array filter. 

GRAPHICAL REPRESENTATION OF TARGETS 

The first task was to represent graphically the amplitude transmittances of a color 
bar chart (red, green, blue, cyan, yellow, and magenta patches), a gray scale, and 
a resolution target. Actually the density levels of all the subsample elements had to 
be specified. This was accomplished directly by visualizing the equivalent encoding 
filter geometry as shown m Figure 40 and writii^ down the resulting densities as in 
Figure 41 . By letting a black subsample patch fall in the upper left-hand comer of 
the blue patch (i. e. , a patch exposed only by blue light) it becomes evident that a 
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vertical ruling pattern must be formed since 
the black and yellow subsamples permit no ex- 
posure and the cyan and magenta subsamples 
transmit blue light that produces moderate 
density levels. Similarly, the red light patch 
consists of a checkerboard pattern since only 
the yellow and noagenta subsamples produce 
moderate densities. Since the yellow and cyan 
subsamples are transparent in green light, the 
green patch consists of a horizontal ruling pat- 
tern. Finally the subtractive colors such as 
yellow have one subsample element of minimum 
density because they transmit two of the three 
primary color components; tiiat is, the yellow 
subsample produces the highest exposure level 
in yellow light. Only three density levels were 
required to specify the color bar chart: maxi- 
mum density, which was set at 3 . 08 because 
higher e:jq)0sure levels would have caused severe distortion (blooming) of the square 
sub sample; moderate density of 0. 62 corresponding to about 50 percent amplitude 
transmittance; and low density of 0.0 corresponding to 100 percent amplitude 
transmittance . 


B - BLACK 
M - MAGENTA 


C - CYAN 
Y - YELLOW 

7 1 06063 


Figure 40. Equivalent Encoding 
Filter Geometry 


^ecification of the gray scale densities as shown in Figure 42 follows the same 
approach except that in this case only neutral exposures of varying intensity increments 
(quarter stops) were simulated. This target required seven different density levels. 

Finally the resolution target densities shown in Figure 43 were determined on the 
basis of a high-contrast black-and-white exposure with the highest frequency 
(4 ^p/mm) equal to one-half of the square array frequency (8 ^p/mm). Only maximum 
and minimum densities were required for this target. 
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Figure 41. Color Bar Chart Densities 
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Figure 42 . Gray Scale Densities 
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PHOTOGRAPHIC FACSIMILE GENERATOR 


The next step in this demonstration was to translate the density geometries of 
Figures 41, 42, and 43 into a set of commands in the form of paper tapes to drive the 
photographic facsimile generator. Normally this step-and-repeat device (Figure 44) 
makes precise exposures of preselected apertures on a photographic plate to generate 
a spatial filter in response to commands from a computer. These commands are 
generated by the computer software program from digital data supplied to the computer. 



OPTICAL-MECHANICAL 
ASSEMBLY COVER 

STAGE COVER 
CONTROL PANEL 


CABINET 


ACCESS DOOR 


LATCH 


6908360 


Figure 44. Photographic Facsimile Generator 
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In this study, however, the device was programmed to print out square array encoded 
color signals. 

A paper tape reader interface unit and a computer buffer interface imit provide 
appropriate signal interfaces between the Photographic Facsimile Generator and the 
signal source. The paper tape reader interface permits the preparation of a complete 
program for reproducing a spatial filter with unique characteristics by separate periph- 
eral equipment and then storing these data on paper tape for future or repeated use . 
Similarly, the conq>uter buffer interface permits direct control of the Photographic 
Facsimile Generator with a Sigma 7 Computer. The conqjuter buffer interface unit 
operates on a priority interrupt with the computer software arranged for accepting 
data from magnetic tape in a Fortran format. Typical specifications for the Photo- 
graphic Facsimile Generator are given below. 

Photosensitive Plate 


Size 

Thickness (glass) 


4 by 5 in. 

1/8 to 1/4 in. 


Step Increments 


Basic unit 
Preset step 

Internal step set range 


0.00025 in. 

0.0025, 0.0050, 0.0100 in. 

Any multiple of 0.00025 times 1 through 127 


Slewing Rate 


0.1 in. /sec 


Aperture Holder 


Accepts 1 in. square material up to 1/4 in. thick. 


Aperture Supplied 


0. 002 in. copper-beryllium photoetched; hole 
corresponds to preset steps and to 0. 020 for 
internal set. 


Magnification 


10. 5X 


Stage 


Type 

Travel (nominal) 


x-y cross slide 

1 in. X and y directions; limit switches set at 
1.1 in. travel 


Image Location 


Centrally located on the 4 by 5 in. glass plate 
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Exposure Duration 

Density 

Lamp 

Lair^ Current 

Operational Modes 

Repetition Cycle 

Exposure Repeatability 

Exposure Uniformity 

Position Accuracy 

Power Requirements 

Voltage 

Current 


20 msec 

Two independent ten-position discs 
Quartz -halogen, 100 W 
7 to 8.5 A 
Manual -automatic 

Five to seven exposures per second 
±2 percent 

Less than 1 percent for a 0, 0025 image 
Better than 0.0004 in. 

117 Vac, single-phase, 60 Hz 
15 A 


DRIVE TAPE PREPARATION 

Before preparing the drive tapes the Photographic Facsimile Generator had to be 
calibrated in terms of output density versus the two density disc positions shown in 
Figure 45. Each disc has ten different levels and in their cascaded arrangement 
produce one himdred different exposure levels. For each element of the printout 
array, (400 by 400) elements, a uniform intensity square (0. 0025 x 0. 0025 in.) is 
exposed to one of the hundred different levels. The elements are displaced by 
(0.0025 in.) to form a continuous array. 

Only seven density disc positions and one no exposure command were required in 
tape preparation; these are summarized in Table 10. Each target was abbreviated 
to about one-third of the format potential in an attempt to save printout time and 
increase realiability. Thus the targets were only 130 rather than 400 spaces wide. 

It was hoped that all three targets could be printed out side by side on a single photo- 
graphic plate, but this task proved too time consuming and error prone. 
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Figure 45. Optical-Mechanical Assembly of Photographic Facsimile Generator 
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Table 10. Density Disc Calibration 


Output Density 

Density Disc Position 

3.08 

50 

0.79 

57 

0.62 

65 

0.58 

72 

0.42 

89 

0.27 

74 

0.12 

64 

0.00 

Command C 


TARGET PRINTOUT 

The paper drive tapes were then used to command the Photographic Facsimile 
Generator which produced the printouts (Figures 46, 47 and, 48) on High Resolution 
Plates that were processed to a gamma of 2. Half a dozen printout errors are 
evident in Figure 46, but these were judged to be only minor and in no way interfered 
with the principle being demonstrated. A spatial phase error is also visible in just 
about the middle of Figure 47, whereas the resolution target of Figure 48 is virtually 
error free. 

Recall that the primary objective has been to demonstrate the feasibility of printing 
out in completely automatic spatial registration on monochrome film, the three chan- 
nels of information required for full color imagery. The irqjut to the facsimile gen- 
erator is the electronic signal from a punched tape taken as the equivalent to a signal 
from an electronic buffer with the three signals in temporal registration. The devel- 
oped printout of the Photographic Facsimile Generator gives the resolution target 
and the gray scale immediately with results superficially similar to a halftone print. 
The printout for the three channel color appears in a black-and-white image, which 
is inserted into a Total Optical Color (TOC) projector to produce the color image. 
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Figure 46. Printout of Color Bar 
Chart (at approx. 5. 4X) 
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Figure 47. Printout of Gray Scale 
(at approx. 5.45 q 
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Figure 48. 


Printout of Reso- 
lution Target 
(at approx. 5.4X) 
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TARGET RETRIEVAL 


The TOC projector used to retrieve the square array targets was similar in 
principle to the projector shown previously in Figures 22 and 23 except that longer 
focal length optics were required to separate the diffraction orders spatially. The 
optical arrangement is shown in Figure 49. Lenses L^^ and Lg are both 495 mm focal 
length lenses used to image at 1:1. A 450 watt Xenon arc was used to provide white 
light illumination. A color retrieval filter of the geometry shown in Figure 50 was 




COLOR FILTERS 
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R - 25o 
B - 47 
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Figure 50. Color Retrieval 
Filter Geometry 
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aligned with the Fourier transform array shown in Figure 51 and placed directly in 
front of the projector lens L^. The color bar chart was retrieved with the results 
shown in Figure 52, which is a Kodacolor print of a negative film exposed directly 


Figure 51. Fourier Transform Array 
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Figure 52. Retrieved Color Bar Chart 
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in the color image plane. The gray scale reconstruction is shown in Figure 53, and 
the equivalent retrieval densities are given in Table 11 with step 2 normalized to 
zero density. 


Table 11. Equivalent Gray 
Scale Densities 


■ Figure 53. Retrieved 
Gray Scale 


7 t 06076 


Step Density 


1 

2 

3 

4 


0.0 

.18 

.28 

.40 


5 .52 

6 .62 

7 1.00 (or greater) 


Finally the resolution target retrieval is shown in Figure 54, and the spatial fre- 
quency progression is given in Table 12. 



Figure 54. Retrieved 
Resolution Target 
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Table 12. Resolution Target 
Spatial Frequencies 


Group Frequency (gp/mm) 


1 

2 

3 

4 

5 

6 


4.0 
2.67 

2.0 
1.6 
1.33 
1.14 
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IMPLEMENTATION 


As emphasized earlier, this part of the program was to demonstrate the feasibility 
of making use of a high speed scanner-printer to print out monochrome information 
from which three-channel color could be retrieved. Since TOC systems for color 
encoding black-and-white films have generally been concerned with color filter grat- 
ings overlaying the film to be exposed, some conceptual changes are required in 
adapting to electronically recorded images. As it develops, the implementation for 
retrieval of images requires little or no modification to make use oT the electronically 
recorded images . The Photographic Facsimile Generator is used to simulate the 
operation of what would generally be expected to be a higher speed device. In its im- 
plementation the square array method would either require a delay-storage device of 
sufficient storage time to hold one line of print and with a resolution time small enough 
to ensure registration as the stored information is printed out, or a dual writing beam 
device. A third approach described in the program proposal, but not demonstrated 
in this study, is to utilize a scanner with a "wobble" trace that prints out in a hexag- 
onal as opposed to a square array. This approach eliminates the requirements of 
either a delay storage or dual beam printing device, but it introduces the somewhat 
unconventional "wobble " beam trace shown in Figure 55 . 
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Figure 55 , Hexagonal Array Modulation Geometry 
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SECTION 5 


CONCLUSIONS 

In this study we have demonstrated a, multiplexing method capable of storing and 
retrieving twelve images at 10 ip/mm resolution with a brightness dynamic range of 
up to 100:1. Any dimensional distortions introduced by the method are inconsequential 
because they are small compared with the system spread function. To maintain com- 
pletely linear tonal reproduction (i.e., gray scale) over the 100:1 output dynamic 
range, a^i intermediate square root storage step is required before exposing the com- 
posite recording. Crosstalk in the retrieval images was more difficult to suppress 
as the number of channels was increased from seven to twelve. For a seven-channel 
recording, crosstalk was vmdetectable subjectively, but it became noticeable for nine- 
and twelve -channel recordings. 

The storage method demonstrated is practical in the sense that it was implemented 
with existing materials and state-of-the-art techniques. Exposure and processing 
controls required to minimize crosstalk are of a lype that can be automated and do not 
inherently require the judgment of a trained technician. The question of whether this 
multiplexing method is acceptable must await an evaluation of this study since we are 
not in a position to determine whether the image qualily is subjectively acceptable. 
Before the economics of this method can be realistically analyzed, a specific system 
must be defined. 

In addition, we have demonstrated a step-and-repeat method for printing out and 
reconstructing three-zone color imagery with monochrome film such that registration 
accuracy of the zones is not degraded. The results provide a positive confirmation 
of the square array modulation theory as well as the method’s registration precision. 
The resolution and dynamic range demonstrated were determined chiefly by the 
Photographic Facsimile Generator characteristics, which were incidental rather than 
fundamental to this simulation experiment. 
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All of the work statement objectives of this study were accomplished, and for the 
most part the results were those generally anticipated at the outset. One unexpected 
results was the superior performance of Recordak AHU Type 5460 processed in a 
direct reversal mode. The moderately wide output dynamic range and the hightened 
crosstalk sensitivity in regions of low average exposure were other interesting re- 
sults that would have been difficult to predict from the theoretical system model. 

Finally, the results of this study indicate that the primary objective is best 
achieved by limiting the number of input channels to seven. 
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SECTION 6 


RECOMMENDATIONS 

Assuming that the image quality of this storage and retrieval method is acceptable, 
we recommend two areas for future work: 

1. Using seven actual ERTS spectral zonal separation negatives (7-1/2 inch 
square format) with registration marks as furnished by groxmd recon- 
struction printers, prepare overlay multiplex recordings with all channels 
in mutual registration. 

2. Design and fabricate a prototype retrieval projector that has the potential 
for relatively low production cost and can be used to demonstrate retrieval 
image quality. 

To pursue high speed scanner-modulator systems for printing out three-zone ERTS 
imagery would be redimdant at the present time, first because the three color laser 
scanner currently under development for NASA will solve the same problem if suc- 
cessful, and second if registration can be adequately maintained in the seven-channel 
recording recommended above, the problem will once again be solved. 

The two areas recommended for future action could be carried out either sequen- 
tially or concurrently since most of the projector design would be independent of re- 
cording preparation. Although projector design is vital in answering the question of 
system economics,it is not however essential in demonstrating large format retrieval 
and registration since that could be accomplished with a bench mock-up similar to 
the approach taken in this study. 
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APPENDIX A 


DESIGN REQUIREMENTS TO EXTEND OVERLAY MULTIPLE 
IMAGE METHOD TO FULL FORIVLAT RECORDINGS 

In this discussion we will restrict our considerations to a seven-channel system 

t 

utilizing Recordak AHU Type 5460 Microfile film because that combination appeared 
most promising in this study. However, most of the design considerations could be 
extended to higher channel capacity if desired. 

COPY CAMERA 

Three major design elements must be addressed in the copy camera; first, the 
ruling must be capable of encoding 7-1/2 x 7-1/2 inch images; second, the copy 
lens must be essentially distortionless and free of vignetting; and third, a suitable 
pressure platen must be used to ensure uniform and controlled ruling contact. 

Typical ruling specifications are given in Table A-1. Other than a general scaling 
of size, the ruling specifications differ from the ruling used in this study only by 
the antireflection coating requirement. 

According to the manufacturer’s specifications, the 160 mm f/5. 6 Fax-Nikkor 
lens used for this study can accommodate an 11 x 11 in. format (see Figure A-1). 

To minimize the effect of lens aberrations and the possibility of distortions and 
vignetting, however specifying a lens with greater format capability would be 
advisable. The 300 mm f/7 Fax-Nikkor lens (Figure A-2) can accommodate a 
50 percent larger format or about 17 by 17 in. , and therefore we recommend that 
lens or one equivalent to it for the 7-1/2 inch format application. 

The pressure platen can be one of many various designs, but the pneumatic bag 
approach sketched in Figure A-3 seems well adapted to this relatively large format 
application. The platen consists of a smooth rubber membrane of fairly high 
density and thickness, which is sealed by means of a frame or vulcanized to a metal 
back plate. A port is provided in the back plate so that a pressurized gas can be 
alternately used to inflate the "bag" and then evacuate it. When inflated the 
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Table A-1. Encoder Rviling Specification for Imagery 7-1/2 Inches Square 


Substrate 


Ruling Frequency 
Duty Cycle 

Transmittance 
(micro) . 

Transmittance 

(macro) 


Ruling Format 

Ruling Material 

Antireflection 

Coating 


Optical quality waterwhite glass of nominal 0,250 in. 
thickness with a flatness tolerances of no more than 
10 Mnges/in. Diameter should be 12 in. ± 1/64 with 
chamfered edges. 

100 fp/mm (±2 percent). 

50 percent square wave (i.e. , bars and spaces equal 
to 0. 0005 ± 5 percent). 

Bar transmittance should be 1 percent or less. 

^ace transmittance should be 94 percent or greater. 

For sample patches greater than 50 periods (i. e. , 

5 mm width or diameter) average transmittance should 
be 47.5 percent ±2,5 percent, but variations over the 
ruling format should not exceed ±1.0 percent 
transmittance. 

Minimum of 11 in. diameter concentric with substrate. 

Metallic (equivalent to chromium or Inconel), 

On the ruling side, reflection must be less than 0.4 
percent over the wavelength range of 400 to 650 nm and 
not exceed 0. 7 percent from 360 to 700 nm. No anti- 
reflection coat IS required on the substrate surface 
opposite the ruling surface. The antiroflection coating 
must be extremely durable and scratch resistant, as 
would be achieved with titanium dioxide, sapphire, and 
fused silica dielectrics applied to hot substrates with 
an electron beam evaporator. No spit or spatter can 
be tolerated, and the surface must be as smooth as 
the original substrate polish with the exception of the 
ruling bars. 



160mm f/5.6 Fax-Nikkor 


SOOmm £/7 Fax-Nikkor 


Focal length 

160mm 

Max. aperture 

f/5.6 

Construction 

6 elements 
4 groups 

Standard magnification 

1 X 

Range of magnification 

1/3 X - 3 X 

Image area 

420mm at f/5.6 

Picture angle 

72* 

Overall working distance 

640mm 

Correction wavelength range 

350m^l— 700m/u 

Vignetting 

0% at f/8 

Aperture scale 

5.6, 8, 11. 16. 22 

Mounts 

Screw 

(d=62mm p«0.75mm) 
& Adapter plate 
(o.d.=77.5mm) 

Dimensions: max. diameter 

70mm 

max. length 

53mm 

Weight 

320g 

Shutter 

B 





Focal length 

300mm 

Max. aperture 

Ml 

Construction 

6 elements 
4 groups 

Standard magnification 

IX 

Range of magnification 

1/3X -3X 

Image area 

630mm 

Picture angle 

52* 

Overall working distance 

1,200mm 

Correction wavelength range 

350m^-700mjLI 

Vignetting 

0% at f/7 

Aperture scale 

Slit for slip-in aperture 
control plate provided 

Mounts 

Screw 

(d=95mm paimm) & 
Adapter plate 
(o.d.a 120mm) 

Dimensions: max. diameter 

103.5mm 

max. length 

96.5mm 

Weight 

1.700g 

Shutter 

B 
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Figure A-1. Copy Lens Speei- Figure A-2. Copy Lens %)ecifica- 

fications Used for 35 nun tions Recommended for 7-1/2 in. 

Format Study Square Format 
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Figure A-3. Cross Section of Pneumatic Pressure Platen 


rubber membrane presses the film into contact with the ruling, and when evacuated 
the pressure is released. The thickness, density, and texture of the rubber membrane, 
as well as the pressure differentials to be used, would have to be verified experimen- 
tally. This approach has the advantage of providing a pressure pad that conforms to 
any film base irregularities, which would tend to present a more serious problem for 
long width films than for the relatively narrow 35 mm chips used for this study. 

FILM AND PROCESSING 

Although 9 in. wide Type 5460 microfilm is a special-order item, we have confirmed 
that Eastman Kodak will supply it in minimum order quantities (20 rolls of 100 ft 
lengths). For a small number of diqjlicate recordings, say five or less, we recom- 
mend the same basic reversal chemistry used in this study except that it should be 
ad^ted to automatic machine processing. For a large number of duplicate copies, we 
recommend that the original multiplexed recording be processed as a relatively low 
gamma (y < 1) negative on either Minicard film, Scientia 10E50 (by Agfa-Gaevert), or 
Type 5460. Release positives would then be contact-printed onto Minicard film pro- 
cessed to a y' < 2 such that yy' = 2. 
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RETRIEVAL PROJECTOR 


The large format retrieval projector can be of the same general design as shown 
in Eigure 23. The condenser system however has to be increased to about 10. 6 in. 
diameter. Using the 300 mm Fax-Nikkor as a projection lens at IT requires that the 
condenser lenses be about 11.3 in. in focal length, which makes the condensers just 
greater than f/1 and therefore practical to fabricate. Such lenses are not off-the- 
shelf items, but since they need not be high qualily lenses (m fact, they do not even 
have to be achromats), they could be fabricated at reasonable cost. In large quanti- 
ties it would probably be economical to use plastic lenses, possibly of a light-weight 
Fresnel design. 

We have already indicated that the projection lens should be a 300 mm Fax-Nikkor 
or an equivalent. The only other significant component that need be considered is the 
light source. To achieve even modest illumination levels over such a large format 
would require a source^ in the 500 watt class. A compact xenon or mercury arc would 
be preferable. A small filament tungsten halogen lamp would provide a more econom- 
ical source if the projector were not used for extended periods; however, retrieval 
brightness would be lower than that achieved with an arc source . 
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The subtractive modulation approach requires that black (blank) be printed in the 
fourth zone . The subtractive approach simplifies not only the film problem, but also 
the scanner since three channels rather than four are printed. 

Since the additive and subtractive colors are related linearly, a taking system 
utilizing either approach can be accommodated. The conversion merely represents 
another digital processing step prior to generating the film recording tapes. 

It should be pointed out that if all the color information is present, the TOC ap- 
proach permits complete flexibility in color balancing in the viewing device. 

The specified format for the final imagery is a 7 , 5 in. square frame at a resolution 
of 10 P,-p/ mm . Since the imagery is sampled, the sampling frequency must be at least 
twice, and preferably as much as three times, the maximum bandwidth of the imagery. 
The sampling frequency, then, must be at least 20 ^p/mm and preferably as high as 
30 ^p/mm. As shown in Figure 8, each sample consists of 4 subsamples. At the 
resolution discussed, each subsample will have a maximum size of 25 x 25 microns 
and preferably will be as small as 16 x 16 microns. 

hi a step-and-repeat printer such as the Photographic Facsimile Generator, each 
printed spot is equivalent in size and shape to the image of an aperture. The image 
of the aperture is positioned at each point to be printed and the shutter opened long 
enough to yield the desired exposure . 

In the ideal case, the aperture image will have sharp edges, and the positional 
accuracy will be such that there is no overlap or underlap between adjacent exposures. 
More realistically, neither situation will be entirely true and there will be some over- 
lap between adjacent color samples. If the film recording is linear, this overlap has 
no effect except to desaturate the color information. Therefore, if the overlap is kept 
small, the effects will be negligible. 

In -a continuous scanner, the ejiposure is the convolution of the writing spot and a 
rectangle with uddth corresponding to writing time. Thus we have 

X 

E(x) = ^ R(t) S(x - vt) dt , 

0 
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APPENDIX B 


DESIGN REQUIREMENTS FOR A HIGH SPEED SCANNER-MODULATOR 
TO IMPLEMENT THREE-CHANNEL INTERLACE STORAGE 

The major considerations in the design of a printer to generate three-color TOC 
(Total Optical Color) are the required resolution, the format size, the dynamic range 
of the recorded image, and the required printout speed. These parameters in turn 
determine the required spot size and shape, the scanning speed and signal bandwidth, 
and the characteristics of the recording medium. 

Since the TOC imagery is to be recovered in a spatial filtering system, careful 
consideration must be given to gamma control to ensure faithful reproduction with 
minimal crosstalk. 

We assume that prior to recording on film, the data will be stored digitally on a 
mass storage medium such as magnetic tape. We further assume that nothing pre- 
cludes the storage of the various color chaimels on the same tape or on two or three 
separate tapes. Any operations such as changing from additive to subtractive modu- 
lation can be performed prior to generating the film recording tapes. 

Since the spatial filtering operation associated with recovering the color imagery 
requires a transparency that is linear in amplitude, the intensity transmission must 
be a square function. Therefore, to ensure linear filtering, the film must be devel- 
oped to a gamma of -2. The further the deviation from a gamma of -2, the more 
nonlinear will be the filtering operation and the more crosstalk that will ensue. 

To compensate for this data expansion, precompressing the data with a square 
root function might be desirable. It is anticipated that this operation would be done 
prior to exposure. 

As pointed out earlier, the colors can be encoded with either subtractive or addi- 
tive modulation, and theoretically these two approaches are equivalent. The additive 
approach, however, requires that the white information, which is the sum of the 
three colors, be printed in the fourth zone. Since the white channel has a wider 
dynamic range than any of the color chaonels, the constraints on the film would be 
quite severe. 
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where E(x) is the exposure on film, R(t) is 1 during the writing time and 0 elsewhere, 
and S(x) is the function describing the writing spot. 

Even if S(x) is rectangular, the exposure will be trapezoidal with tails equal to one 
half the spot width. Therefore, either overlap or gaps between adjacent samples must 
be tolerated. Roughly speaking, if the overlap is to be less than 20 percent, the half 
spot width must be less than 20 percent of the equivalent writing time. Since, by 
previous calculations, the subsample spot size is between 16 and 25 microns, the 
writing spot must be on the order of 7 to 10 microns. 

The effects of spot size versus color quality are not readily determined theoretically. 
Although the figure given may be conservative, it should represent a lower limit on 
spot size. 

Based on the restriction on spot size, consideration can be given to the printing con- 
cept. With a format of 7-1/2 inches and a subsample size of 16 to 25 microns, the reso- 
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lution elements per line are on the order of 10 . This figure precludes a CRT printer. 

The two most likely approaches, depending on speed requirements, are a step-and- 
repeat printer similar to the Tech/Ops Photographic Facsimile Generator or a rotating 
mirror laser scanner. 

The Tech/Ops Photographic Facsimile Generator operates at up to 20 exposures per 
second, but this type of device could be made to operate at up to 100 ejq) 0 sures per 

second. Even at 100 exposures per second, a 7-1/2 in. square format consists of more 
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than 5 X 10 exposures and would therefore take longer than 150 hours (5x10 sec) to 

print. 

A laser printer operating at an average data rate of 1 MHz would require about 
1 minute to print the same format. With a duty cycle of less than 100 percent, the 
peak data rate would have to be correspondingly higher or the printing time longer. 

The last item to be considered is the dynamic range of the recorded image. A nom- 
inal requirement of 16 density levels over a range of 0. 0 to 2. 0 is well within the capa- 
bilities of film, lasers, and laser modulators. The density range requires an extinc- 
tion ratio of the modulator of 100:1, which is well within the state of the art. Sixteen 
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equal levels correspond to density increments of 0.125, which in turn correspond to 
changes in transmission of about 5.5 percent. The two most serious problems in 
attaining this figure are modulator drift and laser noise. The mode switching noise 
in HeNe lasers can be limited to 1 percent or less. By using feedback techniques, 
both modulator and laser drift can be held to well within 1 percent. 

A printing system that will yield TOC coded color imagery at the rate of about one . 
frame per minute is apparently well within the state of the art. We assume that pre- 
processing can be done digitally and the three-color information stored on magnetic 
tape or a similar medium. The most likely printing technique uses a rotating mirror 
laser printer with a spot size of about 10 microns. A medium performance modulator 
in conjunction with a HeNe laser provides sufficient d 3 mamic range for good quality 
color retrieval. 
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